Binary black hole mergers: large kicks for generic spin orientations 
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We present results from several simulations of equal mass black holes with spin. The spin magni- 
tudes are S/m 2 = 0.8 in all cases, but we vary the spin orientations arbitrarily, in and outside the 
orbital plane. We find that in all but one case the final merged black hole acquires a kick of more 
than 1000 km/s, indicating that kicks of this magnitude are likely to be generic and should be ex- 
pected for mergers with general spin orientations. The maximum kick velocity we find is 2500 km/s 
and occurs for initial spins which are anti-aligned in the initial orbital plane. 

PACS numbers: 04.25.Dm, 04.30.Db, 04.70.Bw, 95.30.Sf, 97.60.Lf 
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I. INTRODUCTION 



Currently several laser-interferometric gravitational 
wave detectors such as TAMA Q, LIGO and GEO [3 
are already operating, while several others are in the 
planning or construction phase One of the most 

promising sources for these detectors are the inspirals and 
mergers of binary black holes. The gravitational waves 
radiated during the inspiral and merger of two black holes 
carry energy, linear momentum and angular momentum 
fluxes. The radiated energy is positive and thus the sys- 
tem will always lose energy. Except in the case of cer- 
tain highly symmetric situations, such as head on colli- 
sions of equal mass non-spinning black holes, the system 
usually also loses angular momentum. Generic systems 
with either unequal masses or spins with no particular 
alignment will also radiate linear momentum, resulting 
in a non-zero residual velocity (also known as recoil or 
"kick") of the final black hole. The magnitude of this 
recoil is important in a variety of astrophysical scenar- 
ios, such as the cosmological evolution of supermassive 
black holes @, i, 0, H ©, El El or the growth and re- 
tention of intermediate-mass black holes in dense stellar 
clusters QJ, EI EE EE EE El El ■ For a binary in almost 
circular orbit, the direction of the instantaneous linear 
momentum flux rotates in the orbital plane with the an- 
gular velocity of the system. Thus, when the binary goes 
through one orbital period, the average linear momen- 
tum flux will be close to zero. The only net effect comes 
from the fact that the inspiral orbits are not perfect cir- 
cles. Most of the kick is thus accumulated during the last 
orbit and subsequent plunge of the two holes, when the 
motion of the two holes is no longer quasi-circular and the 
averaged linear momentum flux is much larger. Several 
analytical estimates of the kick velocity have been pub- 
lished in recent years [IE HH, HE [23| . All these estimates 
have been derived using approximations (such as post- 
Newtonian theory) which break down during the last or- 
bit and the subsequent merger. Hence these calculations 
contain a high degree of uncertainty. In addition, there 
are calculations which employ the close limit approxima- 
tion [2j|, |25[. However, these calculations can only model 
well the plunge and miss contributions from the last or- 
bit. Furthermore, results from several numerical simula- 



tions have been published [26|, |27j, \M |29j, |30|, |31|, [3J, |33j 
recently. The resulting maximum kick velocities found 
so far are: 175 km/s for non-spinning unequal mass 
black holes [28|, and 2500 km/s for equal mass black 
holes with anti-aligned spins S/m 2 ps ±0.8 in the or- 
bital plane [3l[ . The latter value is so high that the final 
black hole would escape e.g. both from dwarf ellipti- 
cal and spheroidal galaxies (with typical escapes veloci- 
ties of below 300 km/s) and from giant elliptical galax- 
ies (2000 km/s) [5j. Extrapolations from simulations of 
black hole binaries with spins in the orbital plane pre- 
dict velocities as high as 4000 km/s 32]. The question 
thus arises if such high kick velocities are generic for 
spinning black holes, or if the results of 2500 km/s or 
even 4000 km/s arise only if the initial spins are exactly 
anti-aligned and confined to the initial orbital plane. To 
answer this question we have performed several simu- 
lations of black hole binaries on initially approximately 
circular orbits with spins in arbitrary directions. We find 
that except for one case the kick velocities are well above 
1000 km/s. Thus, kicks of this magnitude seem to be a 
generic feature of black hole mergers with spin around or 
above S/m 2 « 0.8 



II. NUMERICAL TECHNIQUES 

Our simulations are performed using the "mo ving 
punctures" method [U [35| with the BAM code [H, [371 
which allows us to use moving nested refinement boxes. 
We use 10 levels of 2:1 refinements. The outer bound- 
aries are located 240M away from the initial center of 
mass (M being the sum of the initial black hole masses), 
and our resolution ranges between 9M on the outermost 
box to M/56.9 near the black holes. Since our simula- 
tions have no particular symmetries we cannot use the 
usual memory saving techniques and simulate only one 
quadrant of the numerical domain. Thus our simulations 
take four times more memory and time than e.g. the ones 
reported in [381 ] . This makes standard convergence tests 
very costly. Thus, the goal of this work is merely to give 
kick velocity estimates which can be used in astrophysics, 
instead of extremely accurate results. For the same rea- 
son, we present a convergence test for only one particular 
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FIG. 1: This plot shows the energy radiated during the merger 
of a binary where both initial spins (S/m 2 = 0.8) are anti- 
aligned to the orbital angular momentum. The upper panel 
shows the energy radiated for the resolutions hi = M/56.9, 
hz — M/61.6, hs = M/66.4. In the lower panel we show the 
differences scaled for fourth order convergence. 



case in order to establish that our numerical setup gives 
results that are in the convergent regime. In Fig. [T] we 
show the energy radiated during the merger of a binary 
where both initial spins have magnitude S/m 2 = 0.8 and 
are anti-aligned to the orbital angular momentum. The 
lower panel of this plot demonstrates that the resolutions 
hi, hi and h% are in the convergent regime. All other re- 
sults which are given below, have been obtained using the 
lowest resolution hi — M/56.9. From this study it is also 
possible to obtain error bars due to the finite differencing 
scheme used here. These error bars are of the order of a 
few percent. Notice however, that there are other sources 
of errors, the biggest of which is the extraction of the ra- 
diated energy and momentum at a finite radius. This 
radius should ideally be far from the sources. However, 
the outer regions of our grid are not very well resolved, 
since the highly refined boxes do not extend far from 
each black hole. We find that the best compromise for 
our grid structure is to extract at r ex ~ 30M. Figure 
[5] shows the radiated momentum in the z-direction ex- 
tracted at different radii. In this example both initial 
spins are perpendicular to each other and lie in the ini- 
tial orbital plane. We use the difference between these 
results at different radii as error bars for the kick veloci- 
ties presented below. 



III. RESULTS 

In order to gauge possible kick sizes we have per- 
formed several simulations with different spin orienta- 
tions. In each case the spin magnitude of each black hole 
is S/m 2 = 0.8 and the initial orbital angular momentum 
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FIG. 2: The z-component of the radiated momentum for the 
-1,135/-1,-135 run. The initial spins are perpendicular to each 
other and lie in the initial orbital plane. This plot shows that 
on our relatively small grids the location of the extraction 
radius plays a big role and is the main contributor to the 
error of ±20% quoted here. 



is along the z-axis. The magnitude of the orbital angular 
momentum is chosen such that the orbits are circular ac- 
cording to post-Newtonian approximations [38, 39]. The 
bare mass parameter m^ of each puncture is chosen such 
that the ADM mass m measured at each puncture is ap- 
proximately 0.5M. Because of these approximations, it 
is clear that the orbits will have small ellipticities. How- 
ever, since we are just interested in obtaining a sample 
of possible kick velocities for different initial conditions 
such imperfections are not critical. 

We label each simulation using seven numbers: [q, D, 
Sx/m 2 , 90° - 0i, (pi, S 2 /m 2 , 90° - 6 2 , fa]. Defining the 
larger black hole as #2, we set q = m^/mi > 1. D and 
Sa/to 2 represent the holes' coordinate separation and the 
spin magnitude of hole A {A = 1, 2). The angles 9 a and 
4>a (given in degrees) determine the direction of each spin 
by their polar angles with respect to a coordinate system 
(x, y, z), where z is the unit vector in the direction of the 
binary's initial orbital angular momentum, x is the unit 
vector in the direction of the initial linear momentum of 
the larger hole and y = z A x [42j]. Since our runs all 
have equal masses and spin magnitudes, they are all of 
the form [1,6M,0.8,90° - 01,^1,0.8,90° - 2 ,fa]. The runs 
in Table [J are thus labeled by the angles alone. 

From Table U we see that almost all initial spin ori- 
entations lead to large kick velocities. The only excep- 
tion is when both initial black hole spins point in exactly 
the same direction (see run -31,0/-31,0). In this spe- 
cial case no net linear momentum is radiated (in agree- 
ment with post-Newtonian predictions [3!|). We find 
that the largest kick occurs for anti-aligned spins that are 
perpendicular to the initial angular momentum J^ DM . 
However, we also find substantial kicks for other ini- 
tial spin orientations, and in each case the kick is well 
above 500 km/s. This means that kicks of this magni- 
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FIG. 3: This plot shows the tracks of the black hole centers 
for the -1,180/89,180 run. Here one of the initial spins is in 
the initial orbital plane, and the other one is perpendicular 
to the orbital plane. One can see that the orbits do not stay 
in one plane. 

tude should be expected for generic black hole mergers 
with spins of magnitude of S/m 2 = 0.8 or above. 

This result can be explained in the following way: it 
is known that the black hole spins significantly influence 
the dynamics of the two black holes. For example, the 
merger can be delayed if both spins are aligned with the 



orbital angular momentum [40j. Or the orbits can be 
non-planar as in the case of a test mass in orbit around 
a spinning Kerr black hole [H| . Since the linear momen- 
tum radiated depends on the orbit shape, the kick veloc- 
ity also depends strongly on spin. In fact, if the spins 
are not parallel to the orbital angular momentum, the 
spins and the orbital angular momentum both precess so 
that the resulting orbits do not stay in one plane and be- 
come much more complex than the familiar Newtonian 
circular orbits. Since the instantaneous linear momen- 
tum flux changes direction in a less regular way than for 
circular orbits, one expects that the averaged linear mo- 
mentum flux is significantly larger than for quasi-circular 
orbits. As an example of such a more complicated orbit 
see Fig. [3] which shows the orbits of both black hole cen- 
ters for the case where one of the spins starts out in the 
initial orbital plane and the other one perpendicular to 
this plane. We can clearly see that the orbits do not stay 
in one plane, and that the kick moves the final merged 
black hole downward out of the initial orbital plane. 

From Figs. [1] and [2] we can see that the bulk of the 
radiated energy and momentum is emitted within a time 
of about AQM. This is also the time during which the 
maximum amplitude of the gravitational waves generated 
by the merger pass through the detector. I.e. almost the 
entire kick is accumulated during the final merger phase. 
Table H] gives estimates for the mass and spin of the final 
black hole. 



IV. DISCUSSION 



We have performed simulations of equal mass black 
holes with initial spins of mag nitude of S/m 2 = 0.8. We 
have tested several initial spin orientations and find that 
all but one leads to kick velocities well above 1000 km/s. 
Thus kicks between 1000 km/s and 2000 km/s seem to be 
a generic feature for spins of magnitude 0.8 with general 
orientations. Higher spins will lead to even larger kicks. 
This means that the final black hole will likely escape 
from dwarf elliptical and spheroidal galaxies (with es- 
cape velocities below 300 km/s), and possibly even from 
giant elliptical galaxies (with escape velocities around 
2000 km/s) Q. This result also has important impli- 
cations for a variety of astrophysical scenarios, such as 
the cosmological evolution of supermassive black holes 
or the growth and retention of intermediate-mass black 
holes in dense stellar clusters. Note, however, that a more 
extensive parameter search should be carried out before 
drawing any definitive conclusions. Furthermore, notice 
that all our simulations are for the equal mass case. It is 
clear that the kick velocity will depend on the mass ratio 
(for two proposed models see [lo, HI]), and we expect 
smaller kicks for the unequal mass case. The probability 



for the large kicks obtained here might be quite low, if 
kick velocities drop sufficiently fast for mass ratios differ- 
ent from one. 

Almost the entire kick is accumulated during the final 
merger phase, meaning that small changes in the initial 
conditions are likely to lead to quite different kick direc- 
tions and magnitudes. The reason is that a small change 
in the initial velocities or separation (and also accumu- 
lated numerical errors in the orbital phase), will signifi- 
cantly alter the merger time (3?| , and thus alter the spin 
orientations just before merger. We therefore expect dif- 
ferent trajectories just before the merger which then will 
lead to different final kicks. Note, however, that this sen- 
sitivity to initial conditions does not change our main 
conclusion that generic spin orientations can lead to very 
large kicks in the equal mass case. 
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Run 


0,90/0,-90 


-1,135/-1,-135 


45,0/-45,180 


44,180/-46,180 


-31,0/-31,0 


-1,180/89,180 


-1,150/89,151 


-1,120/89,119 


M^ UM /M 
J* DM /M 2 


0.985 


0.985 


0.985 


0.985 


0.985 


0.984 


0.984 


0.984 


0.842 


0.882 


0.839 


0.881 


0.756 


1.016 


1.016 


1.016 


ra/M 


0.50 


0.50 


0.50 


0.50 


0.50 


0.50 


0.50 


0.50 


m b /M 


0.30 


0.30 


0.30 


0.30 


0.30 


0.30 


0.30 


0.30 


D/M 


6.000 


6.000 


6.000 


6.000 


6.000 


6.000 


6.000 


6.000 


P/M 
L/M 2 
Si x /m 


0.140 


0.140 


0.140 


0.140 


0.147 


0.133 


0.133 


0.133 


0.842 


0.840 


0.839 


0.839 


0.881 


0.798 


0.798 


0.798 





-0.566 


0.566 


-0.576 


0.684 


-0.800 


-0.693 


-0.400 


ci I 2 

Siy/m 


0.800 


0.566 














0.400 


0.693 


S lz /m 2 





-0.010 


0.566 


0.556 


-0.414 


-0.010 


-0.010 


-0.010 


n 1 2 

S 2 x/m 





-0.566 


-0.566 


-0.556 


0.684 


-0.010 


-0.009 


-0.005 


S2 V /m 2 


-0.800 


-0.566 














0.005 


0.009 


S2z/m 2 





-0.010 


-0.566 


-0.576 


-0.414 


0.800 


0.800 


0.800 


M f /M 

JjM 


0.95 


0.95 


0.95 


0.952 


0.96 


0.94 


0.94 


0.94 


0.67 


0.72 


0.68 


0.73 


0.64 


0.81 


0.80 


0.80 


v (km/s) 


2500 


1350 


580 


1350 





1850 


2100 


1950 


±500 


±250 


±150 


±250 


±50 


±350 


±450 


±400 



TABLE I: Initial parameters and resulting kick velocities. The runs are labeled by the spin orientation angles, which for each 
hole are the angle between the spin and the initial orbital plane and the angle between the projection of the spin into the 
orbital plane and the initial linear momentum of hole 2. Here rrib is the bare mass parameter of each puncture and m the 
ADM mass measured at each puncture. The black holes have coordinate separation D, with puncture locations (0, ±D/2,0) 
and linear momenta (=pP, 0,0), so that the initial orbital angular momentum is in the z-direction. The spins have a magnitude 
of S/m 2 = 0.8 and point in the various directions listed here, so that the total angular momentum is not necessarily along the 
2-direction. We also list the estimates for the initial ADM mass M^, DM and ADM angular momentum J^ DM . Furthermore, 
we give the final black hole mass M/ and angular momentum J/ , as well as the kick velocities v for the different cases. 



National Computational Science Alliance under Grants performed at the Charles E. Schmidt College of Science 
PHY060021P and PHY060040T. Simulations were also computer cluster Boca 5. 



[1] TAMA, http://tamago.mtk.nao.ac.jp/. 

[2] LIGO, http://www.ligo.caltech.edu/. 

[3] GEO, http://www.geo600.uni-hannover.de/. 

[4] B. Schutz, Class. Quantum Grav. 16, A131 (1999). 

[5] D. Merritt, M. Milosavljevic, M. Favata, S. A. Hughes, 

and D. E. Holz, Astrophys. J. 607, L9 (2004), astro- 

ph/0402057. 

[6] M. Boylan-Kolchin, C.-P. Ma, and E. Quataert, Astro- 
phys. J. 613, L37 (2004), astro-ph/0407488. 

[7] Z. Haiman, Astrophys. J. 613, 36 (2004), astro- 
ph/0404196. 

[8] P. Madau and E. Quataert, Astrophys. J. 606, L17 

(2004), astro-ph/0403295. 
[9] J. Yoo and J. Miralda-Escude, Astrophys. J. 614, L25 
(2004), astro-ph/0406217. 
[10] M. Volonteri and R. Perna, Mon. Not. Roy. Astron. Soc. 

358, 913 (2005), astro-ph/0501345. 
[11] N. I. Libeskind, S. Cole, C. S. Frenk, and J. C. Helly, 
Mon. Not. Roy. Astron. Soc. 368, 1381 (2006), astro- 
ph/0512073. 

[12] M. Micic, T. Abel, and S. Sigurdsson (2005), astro- 
ph/0512123. 

[13] M. C. Miller and D. P. Hamilton, MNRAS 330, 232 
(2002). 

[14] M. C. Miller and D. P. Hamilton, Astrophysical Journal 

576, 894 (2002). 
[15] H. Mouri and Y. Taniguchi, Astrophysical Journal 566, 

L17 (2002). 



[16] H. Mouri and Y. Taniguchi, Astrophysical Journal 580, 
844 (2002). 

[17] K. Gultekin, M. C. Miller, and D. P. Hamilton, Astro- 
phys. J. 616, 221 (2004), astro-ph/0402532. 

[18] K. Gultekin, M. Coleman Miller, and D. P. Hamilton, 
Astrophys. J. 640, 156 (2006), astro-ph/0509885. 

[19] R. M. O'Leary, F. A. Rasio, J. M. Fregeau, N. Ivanova, 
and R. O'Shaughnessy, Astrophys. J. 637, 937 (2006), 
astro-ph/0508224. 

[20] A. G. Wiseman, Phys. Rev. D 46, 1517 (1992). 

[21] M. Favata, S. A. Hughes, and D. E. Holz, Astrophys. J. 
607, L5 (2004), astro-ph/0402056. 

[22] L. Blanchet, M. S. S. Qusailah, and C. M. Will, Astro- 
phys. J. 635, 508 (2005), astro-ph/0507692. 

[23] T. Damour and A. Gopakumar, Phys. Rev. D73, 124006 
(2006), gr-qc/0602117. 

[24] C. F. Sopucrta, N. Yunes, and P. Laguna, Phys. Rev. 
D74, 124010 (2006), astro-ph/0608600. 

[25] C. F. Sopuerta, N. Yunes, and P. Laguna, Astrophys. J. 
656, L9 (2007), astro-ph/0611110. 

[26] F. Herrmann, D. Shoemaker, and P. Laguna (2006), gr- 
qc/0601026. 

[27] J. G. Baker et al., Astrophys. J. 653, L93 (2006), astro- 
ph/0603204. 

[28] J. A. Gonzalez, U. Sperhake, B. Bruegmann, M. Han- 

nam, and S. Husa (2006), gr-qc/0610154. 
[29] F. Herrmann, I. Hinder, D. Shoemaker, P. Laguna, and 

R. A. Matzner (2007), gr-qc/0701143. 



■5 



[30] M. Campanelli, C. O. Lousto, Y. Zlochower, and D. Mer- 

ritt (2007), gr-qc/0701164. 
[31] J. A. Gonzalez, M. D. Hannam, U. Spcrhake, B. Brug- 

mann, and S. Husa (2007), gr-qc/0702052. 
[32] M. Campanelli, C. O. Lousto, Y. Zlochower, and D. Mer- 

ritt (2007), gr-qc/0702133. 
[33] J. G. Baker et al. (2007), astro- ph/0702390. 
[34] M. Campanelli, C. O. Lousto, P. Marronetti, and 

Y. Zlochower, Phys. Rev. Lett. 96, 111101 (2006), gr- 

qc/05 11048. 

[35] J. G. Baker, J. Centrella, D.-I. Choi, M. Koppitz, and 
J. van Meter, Phys. Rev. Lett. 96, 111102 (2006), gr- 
qc/0511103. 

[36] B. Bruegmann, W. Tichy, and N. Jansen, Phys. Rev. 
Lett. 92, 211101 (2004), gr-qc/0312112. 



[37] B. Bruegmann, J. Gonzalez, M. Hannam, S. Husa, 
U. Sperhake, and W. Tichy (2006), gr-qc/0610128. 

[38] P. Marronetti, W. Tichy, B. Bruegmann, J. Gonza- 
lez, M. Hannam, S. Husa, and U. Sperhake (2007), gr- 
qc/0701123. 

[39] L. E. Kidder, Phys. Rev. D52, 821 (1995), gr- 
qc/9506022. 

[40] M. Campanelli, C. O. Lousto, and Y. Zlochower, Phys. 
Rev. D74, 041501 (2006), gr-qc/0604012. 

[41] M. Campanelli, C. O. Lousto, Y. Zlochower, B. Krishnan, 
and D. Merritt (2006), gr-qc/0612076. 

[42] In this paper, we constructed our initial data sets such 
that the numerical grid coordinate system (x, y, z) coin- 
cides with (x, y, z) 



